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Abstract 

I report on lattice QCD calculations that study the properties of the Oq and 
/o mesons. 

Dh 1 Introduction 

(— I 

I review the recent lattice results for the light J^*^ = scalar mesons. 
^ I The interpretation of many 0+^ mesons in terms of quark and glue degrees 

K>- I of freedom is still not clear [1,2]. The 0"''^ mesons potentially contain glueball, 

tetraquark, meson molecule or even quark-antiquark degrees of freedom. I 
Tij- ■ have recently written a review [3] of light meson spectroscopy from lattice 



O 



(N 



QCD, that contains more detail on many of the topics covered here. 



^ ■ 1.1 Background to lattice QCD 

The physical picture behind lattice QCD calculations is that an interpolating 
^ i operator creates a hadron in the QCD vacuum and after a specific time 

^ I interval the hadron is destroyed. The choice of interpolating operator is 

particularly important for hadrons where it is not clear how the hadron is 
built out of quarks and gluons. 

For example, to create a light flavour singlet 0+"^ hadron, possible inter- 
polating operators are 

01 = qq (1) 

02 = qimQlr>q (2) 

03 = Uj,ia, (3) 
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Mathur et al. [13] 





1.41 ±0.12 


SCALAR [10] 
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Table 1: Lightest strange-light O"*" meson from lattice QCD. 

where Upiaq is a spatial plaquette of gauge fields with O"*""*" symmetry, and q 
is a light quark operator. 

The majority of recent lattice QCD calculations include the dynamics 
of sea quarks and have pion masses as low as 300 MeV [3]. The results I 
will present for scalar mesons largely use the last generation of lattice QCD 
calculations that are quenched or dynamical QCD calculations with pion 
masses above 500 MeV [4,5]. 

There are a number of reasons that lattice calculations of the light scalar 
mesons are challenging. The lattice QCD correlators for scalar mesons are 
more noisy than for p and tt mesons. The light scalar mesons decay via 
S-wave decays, and current lattice QCD calculations are in the quark mass 
regime where some decay channels to two mesons are open. 

Eventually, the issue of dealing with resonances in lattice QCD will be 
dealt with by Liischer's formalism [6] that produces scattering phase shifts. 
This year Liischer's technique for resonances was applied to the p meson for 
the first time, by the CP-PACS collaboration [7]. 

1.2 The flavour non-singlet 0++ and 0+ mesons. 

Although I am going to loop through the lattice results for the lightest O"*""*" 
and O"*" mesons, it is important to classify the states into SU3 multiplets or 
some other classification based on tetraquarks for example. 

In the PDG the lightest strange-light O"*" meson is the i^o(1430) [8]. There 
have also been claims that experimental data is consistent with O'^ / = 1/2 
meson called the k with a mass of 660 MeV [9]. The existence of the k is 
controversial, see [1,2,10] for a discussion. 

In table [1] I collect results for the mass of the lightest 0^ sg meson from 
lattice QCD calculations. The lattice results in table [U are consistent with 
experimental mass of the i^Q(1430), but mostly miss the controversial k par- 
ticle. All the lattice calculations used sq interpolating operators, so may have 
missed the k state, if it is mostly a tetraquark state, with no overlap with sq 
interpolating operators. 

Experimentally the lightest 1 = 1 0++ mesons are the ao(980) and the 
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Burch et al. [15] 





~ 1.45 


Hart et ah [16] 


2P 


1.0(2) 


Prelovsek et ah [11] 
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1.58(34) 


Prelovsek et ah [11] 


2P 


1.51(19) 


Mathur et al. [13] 





1.42(13) 



Table 2: A collection of results from lattice QCD for the mass of the lightest 
non-singlet O"*"^ meson. The P stands for a partially quenched analysis. 

ao(1450) [8]. There have been speculations that the ao(980) meson is a 
molecule or tetraquark state [1,2], so it is interesting to see whether lat- 
tice QCD calculations with qq interpolating operators couple to the ao(980) 
meson. In quenched QCD there is a ghost contribution [14], due to the rjn 
contribution, to the scalar correlator that needs to be subtracted off the lat- 
tice data. I collect together some recent results for the mass of the light O"*"^ 
meson from lattice QCD in table [21 I only include quenched data where the 
TjTT contribution has been corrected for [14]. 

McNeile and Michael [12], in an unquenched lattice QCD calculation 
focused on the mass difference (in the hope that systematics cancel), be- 
tween the l"*"" and the 0"*"+ mesons. The lattice calculation used gauge 
configurations from UKQCD [4] and CP-PACS [5]. The mass of the light 
1+^ state was always higher than the O"*"^ meson. The final result was 
mfej — niag = 221(40) MeV, compared to the experimental result of 245 MeV. 
Lang et al. recently reported masses for the lightest flavour non-singlet O"*"^ 
consistent with the mass of the ao(980) meson, from an unquenched lattice 
QCD calculation using chirally improved fermions [17]. 

The previous lattice QCD calculations, discussed in this section, were in 
a regime where the quark masses were large enough that the decay Oq — > rjn 
was forbidden. Now I discuss the new lattice QCD calculations where the 
decay ao — > rjn is energetically allowed. 

The MILC collaboration [18] originally claimed that they had evidence 
for ao decay to tit] from their calculations with improved staggered fermions. 
Other decays are discussed in [19]. Later work by the MILC [19] and UKQCD [20] 
collaborations showed that the lightest state in the flavour non-singlet 0++ 
channel was actually below the irrj threshold, with improved staggered fermions. 
This was puzzling, because experimentally the ao tttt decay is forbidden 
by G parity. 

In [21], Prelovsek explained the behaviour of the flavour non-singlet O"*"^ 
correlator with improved staggered fermions using staggered chiral pertur- 
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Figure 1: Mass of hghtest state in 0"'""'" channel with 7rri2 decay threshold. 

bation theory. Bernard, DeTar, Fu, and Prelovsek [22] extended the original 
analysis by Prelovsek, and also applied it to the flavour singlet /o meson. 
A larger study, with more sea quark masses, is required to say something 
specific about the mass of the lightest oq meson. 

The ETM collaboration have preliminary results for the mass of the light 
0"'"'^ meson from a nf=2 unquenched lattice QCD calculation with twisted 
mass fermions [23-25] . In figure [1] I plot the mass of the light O^"*" meson 
and the vr + 772 decay threshold as a function of the square of the pion mass. 
There was a bug in the original preliminary analysis of the ao masses from 
ETMC, however the plot in figure [T] is from larXiv:0906.4720l and is correct. 
The mass of the 7]2 was computed by Michael and Urbach [24] . Figure [T] 
shows some evidence for the mass of the 0"^"*" tracking the vr + 772 threshold, 
or at least for it being an open decay channel. 

Some caution is required in the interpratation of the results, we are only 
just starting to deal with mesons with open decay channels in unquenched 
lattice QCD calculations. There is a of order 250 MeV difference between 
the mass of the lightest flavour singlet pseudoscalar meson in lattice QCD 
calculations with n/ = 2 and nj = 2 + 1 sea quark flavours [3] and this will 
be important for the decay thresholds. 

There are other quantities, other than masses, that can help determine 
the quark and glue content of scalar mesons. For example, Narison [26] 
proposed to use the leptonic decay constant of the non-singlet 0"'""'" mesons to 
determine the structure of the Oq meson. The fa^ decay constant of the light 
flavour non-singlet O"*"^ meson has been computed using unquenched lattice 
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QCD [12]. 



(0 I qq\ao) = Majao 



(4) 



See [12,27,28] for a further discussion of this decay constant and the connec- 
tion with the electroweak current. 

A molecule of two mesons should have a very small "wave- function" at 
the origin, hence fag should be small. The definition of is similar to that 
of the pion decay constant. Hence we mean "small" relative to 130 MeV. 
The other measured decay constants of pseudoscalar mesons are within a 
factor of 2.5 to the pion decay constant [8]. The only exception is the decay 
constant of the 7r(1300) that is suppressed [29,30]. A large value for decay 
constant fag does not rule out a qqqq multi-quark meson. 

Using gauge configurations from UKQCD and CP-PACS, McNcilc and 
Michael computed fa^ ~ 480 MeV. Sum rule and model estimates find fag 
in the range 290 to 440 MeV [26,28,31,32]. The fa^ decay constant depends 
on the scale and this should be specified for a more detailed comparison. 

Computing the decay width of a hadron is also very a valuable way of 
identifying a state on the lattice. In [12], it was reported that the experimen- 
tal hadron coupling for the decays ao(980) KK and ao(1450) KK were 
0.9 and 0.5 respectively. A lattice calculation [12] found that the lightest 
hadron in the 0+''" correlator had a coupling to KK of ~ 1, thus providing 
additional evidence that the hghtest state was the ao(980). 

Pennington [33] has recently extracted the two photon decay width of the 
cr from experiment to be r(a 77) ~ 4 kcV. Pennington notes that value of 
r(cr — > 77) can depend quite sensitively on the quark content of the cr [33]. 
Recently a formalism to compute two photon widths on the lattice has been 
developed [34]. Dudek and Edwards have computed r(xo — ^ 77) = 2.4 ± 1.0 
keV, from a quenched QCD calculation [34]. It would be interesting to do a 
similar calculation for light scalar mesons. 

1.3 Flavour singlet 0++ mesons 

The spectrum of the light flavor singlet 0++ mesons is where the 0++ glueball 

is thought to be hiding out. The lightest flavor singlet O"*"^ mesons listed in 
the PDG [8] are: /o(600), /o(980), /o(1370), /o(1500), and /o(1710). There 
are claims that the /o(980) is a molecule or tetraquark [8], so it may not 
couple to qq interpolating operators. 

Morningstar and Peardon [35] obtained Mq++ = 1730(50) (80) MeV for 
the mass of the lightest O"*"^ glueball from quenched QCD. Chen et al. [36] 
recently found Mq++ = 1710(50) (80) MeV. The quark model predicts that 
there should only be two 0+"*" mesons between 1300 and 1800 MeV, so if the 
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mixing between the glueball and qq operators is weak, then the 0"'"+ glueball 
is hidden inside the /o(1370), /o(1500) and /o(1710) mesons. 

Weingarten and Lee [37] used quenched lattice QCD to estimate the mix- 
ing matrix between the glue and qq states. Weingarten and Lee [37] predicted 
that the /o(1710) meson was 74(10)% 0+^ glueball, and hence the mixing be- 
tween the 0"*"+ glueball and qq states was weak. 

There are claims [38] that continuum phenomenology is more consistent 
with a a sizable contributions from the 0"^"^ glueball to the /o(600) and 
/o(980) mesons. 

The SESAM collaboration studied the glueball spectrum on unquenched 
lattices [39]. McNeile and Michael studied the light 0+"'' spectrum with un- 
quenched QCD [40] at a coarse lattice spacing and found the mass of the 
lightest flavour singlet 0++ meson was very light. Using 0"^+ glueball opera- 
tors. Hart and Teper [41] found that 



at a fixed lattice spacing of 0.1 fm. The UKQCD collaboration [20] separately 
studied 0+^ glueball and 0++ qq operators on improved staggered gauge 
configurations, however higher statistics and an analysis similar to the one 
by Bernard et al. is required [22]. 

Unfortunately, the existing unquenched lattice QCD calculations of the 
flavour singlet 0++ mesons don't have the range of lattice spacings where a 
continuum extrapolation can be attempted. In quenched QCD it was found 
that the lattice spacing dependence of the mass of the O"*"^ glueball was 
strong. The use of a Symanzik improved gauge action by Chen et al. [36] 
and, Morningstar and Peardon [35], produced a slightly smaller slope with 
lattice spacing of the scalar O"*""*" glueball mass, than for calculations that used 
the Wilson plaquette action. This is relevant to unquenched calculations, 
because any suppression of the mass of the flavour singlet O"*""*" mass may be 
due to lattice spacing effects. 

The SCALAR collaboration [42], used unquenched lattice QCD, with 
Wilson fermions and the Wilson gauge action, to study the 0++ mesons. At 
a single lattice spacing a ~ 0.2 fm, with qq interpolating operators only, they 
obtain rriqq ~ rUp. The lattice spacing dependence of this result needs to 
quantified. 

In unquenched QCD, both glue and qq states will couple to singlet 0++ 
mesons, so it is better to do a variational fit with both types of operators 
as basis interpolating operators. The variational technique analysis of the 
singlet O"*"^ mesons was done by Hart et al. [43]. A combined fit to 0^+ glue 
and qq interpolating operators with two types of spatial smearing sources 
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Figure 2: Summary of unquenched results for hghtest flavour singlet 0"^"*" 
mesons from [43]. The unquenched results are from SESAM [39], UKQCD- 
I [40], and UKQCD-II [43]. 

was done. The calculation used the non-perturbative improved clover action 
at a single lattice spacing [4]. Configurations from CP-PACS [5] with the 
Iwasaki gauge action and tadpole improved clover action were also used in the 
analysis, because this calculation should be less affected by lattice artifacts. 
A summary plot of the results, in units of ro (I/tq ~ 400 MeV) is in figure [2] 
(updated from [43]). The data with the bursts and squares (with the pion 
masses written near them) in figure [2] shows an additional reduction of the 
mass of the O^"*" state over the pure glueball operators, as used by Hart and 
Teper [41]. 

Mathur at al. [13] recently claimed to get a result for the mass of the 
/o(600) (a) from quenched lattice QCD with pion masses as low as 180 MeV. 
Using the interpolating operator ipj^ipipj^ip they obtain mj(,(6oo) ~ 550 MeV. 
The key part of this work is a three state fit {tt{p = 0)7i{p = 0) , /o(600), 
7r(p = ^)n{p = ^^) using the Bayes adaptive curve fitting algorithm [44]. 
They studied the finite volume effects to distinguish the signal for the res- 
onance from the tttt scattering states [45]. Mathur et al.'s [45] calculation 
is discussed in shghtly more detail in [3]. The effect of sea quarks on this 
calculation needs to be quantified. 

There has also been a recent quenched QCD study [46] of light 0"*"+ states 
with qqTjq interpolating operators that did not see resonant states in the quark 
mass regime they explored. 
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In [43] an attempt was made to compute the decay width for /o decay to 
two pions. Unfortunately much higher statistics will be required to obtain 
an accurate value for that width. 

2 Conclusions 

There is still no consensus as to whether qq operators in lattice QCD calcu- 
lations are coupling to the ao(980) meson. To clear up the many questions 
about the spectrum of the 0+"*" scalar mesons, unquenched lattice QCD calcu- 
lations with tetraquark interpolating operators are required. There is "some" 
evidence that the flavour singlet O"*""^ interpolating operators, in unquenched 
lattice QCD calculations, are coupling to states around or below 1 GeV [43]. 
Although a continuum extrapolation is required for definite results. Recent 
lattice QCD calculations that include the dynamics of the sea quarks are 
working with light enough quarks that the two meson decays of some scalar 
mesons are allowed [3,24,25]. 
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